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Fig. 1-Illustration ol geometry and physics of flow 
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TWO-EQUATION TURBULENCE MODELS 


Ml = q, p k'^/l 


°b G r = -It <p!t £> + p< c *- £ > 


dxj Trfc c)x 


T = si C +r ' Cl c p0k - f2 C2 p r + Cjp ir 

whcr, 


St an dard ^ q 3 Ci= 1.44, C2=1.92, C3=0.0, fi=1.0. 12=1.0, and V=1 .0 

^tended k e^modd. ^ Cl=1 15 c 2 =1.9, C3=0.25 . fi=1.0. f 2 =i.O, and f p =1 .0 

Lam-Bremhorst low-Re model: 2 

Prk = 1.0, Pr e =1.3, Ci=1.44, C 2 =1.92, C3=0.0, fi=(1+0.05/f M ) 3 , t 2 =1-e R t. 
and ^=(1 -e 0 01 65 R k) 2 (1+20.5/R,), where R k =k 1/2 y p/|i and R t =k 2 p/p e 
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Fig. 2(c)-Elfect of y+ on the cross-stream velocity near leading wall 
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Fig. 4(b)-Comparlson of model results with data on trailing wall 



Fig. 4{c)-Comparlson of model results with data on leading wall 
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Fig. 6(a)*Plots o< Nu-ratlo and k-ratlo from the two models 
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Fig. 7(a)<Comparlson on both leading and trailing walls 
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Fig. 8(a}-Comparlson of Eke and low-Re results on leading wall 
(Re=1 0000,Ro=.088,Case B) 



Fig. 8(b)-Comparlson of Eke and low-Re results on leading wall 
(Re=500Q,Ro=.l76,CaseB) 
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Fig. 8(c)-Comparison ol Eke and low-Re results on leading wail 
(Re=5000,Ro=.176,CaseC) 
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Fig- 9-Velocity vectors in the yz-plane 




Table 1 Prediction of Flow Reversal Near the Leading Walt 


Ro 

AT/Tw 

RA1 

Re 

Gr/Re 2 

Row Reversal ? 

0 12 

0.07 

49 

25000 

0.05 

No 


0 13 



0.09 

No 


0 23 



0.16 

No 


0 36 



0.26 

No 


0.48 



0.34 

Yes 


0.07 

196 


0.20 

No 


300 


0.30 

Yes 

0 24 

0.07 

33 

25000 

0.13 

No 


49 


0.20 

No 



196 


0.77 

Yes 



31)0 


1.18 

Yes 


0.13 

49 


0.36 

Yes 


0 16 



0.45 

Yes 


0.23 



0.65 

Yes 


0.07 


12500 

0.20 

No 


0 13 



0.36 

Yes 


0.23 



0.65 

Yes 

0.34 

0.13 

0.16 

49 

25000 

0.73 

0.91 

Yes 

Yes 


0.23 



1.30 

Yes 

0.48 

0.13 

49 

25000 

1.45 

Yes 


SUMMARY 


1 Near-wall grid size has a significant effect on the heat transfer calculations 
when the "wall function" treatment is used. Numerical experiment on the data 
of Morris et al. (1991) suggests that a y+ value in the range of 12 to 42 or so 
yields more accurate results. 

2 The extended k-e turbulence model, while yielding heat transfer results virtually 
the same as those of standard k-e model for low rotation-number flows, 
provides an improvement over the standard k-e model by up to 15% or so in 
heat transfer predictions for high rotation number flows. 

3. Wall-function k-e models predict lower (than data) heat transfer at the trailing 
wall and higher at the leading wall. The need to properly represent the effect of 
rotation in the k-e model equations is realized. 

4 The low-Reynolds number model utilizes a large number of cells and the 
convergence rate is very slow in comparison to the high-Reynolds number 
model using wall function. It is difficult and expensive to obtain a well 
converged solution with the low-Re turbulence model. 
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5. The poor agreement of the low-Re model results with the data makes the low- 
Re mode as an unattractive choice for heat transfer computations in rotating 

number (25000) ° W 31 h ' 9h Rotation number (> 0.24) and high-Reynolds 


6. The extended version of high-Reynolds number turbulence model in 
conjunction with wall function yields satisfactory results for flows with 
isothermal walls as well as uneven wall temperatures. The agreement is within 
5-25% of the data with uneven wall temperatures for flows at Reynolds 
numbers 10000 or higher. 


7. For flows at Reynolds number 5000 or lower, the low-Re model predictions are 
better, especially for the case of uneven wall temperature conditions. 

8. The centrifugal buoyancy may cause a flow reversal near the leading wall 
depending upon the geometry and flow parameters such as rotation number, 
temperature ratio, mean radius ratio and Reynolds number. For the sguare- 
section channel considered here, a criterion of Bo=Gr/Re2 higher than 0 3 is 
predicted to cause flow reversal near the leading wall for flows at Reynolds 
number up to 25000. 
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